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An evaluation of charged pion polarizability using correlation functions measured on a 16 3 x 24 
lattice in the context of the Das, Mathur, Okubo (DMO) sum rule is carried out. The calculation 
is limited to the so-called intrinsic part of the polarizability. This contribution, equivalent to a Eu- 
oo : clidean time integral over vector and axial vector momentum-differentiated propagators, is evaluated 

in both a continuous and discrete sense. In the continuous case, the time behavior of the correlation 
functions is fit to the continuum quark model; the discrete case is handled by a straightforward 
application of Simpson's rule for integration after subtracting the pion contributions. A comparison 
of the implied vector meson and pion decay constants with phenomenological values is carried out. 
Results for the intrinsic polarizability are extrapolated across four quark mass values to the chiral 
■ limit. An extensive discussion of the lattice systematics in this calculation is given. 
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I. Introduction 



Charged pion polarizability is a fascinating laboratory for strong interaction physics. It sits at the crossroads of 
experiment, dispersion relations, sum rules and chiral model results. It also can be evaluated on the lattice and 
turns out to be a sensitive barometer of such inputs as lattice scale, ground state mass values, and models for fitting 
propagator data. It is therefore both an excellent testbed for investigating gauge lattice systematics as well as an 
extremely interesting and fundamental dynamical quantity in itself. 

Some previous experimental and theoretical results on charged pion polarizability (a^±) are presented in Table 
I. Both experimental and theoretical disagreements arc evident. It is thus an opportune time to initiate lattice 
studies which attempt to calculate this quantity from first principles p|]lC[ . Excellent reviews of the experimental and 
theoretical situations are available fTl]] . 

External field methods have been used previously to measure neutral particle polarizabilities in the context of 
lattice QCD | jl2| . Of course these methods would be very difficult or impossible to use in the case of charged hadrons 
because charged particles accelerate in an electric field. The evaluation of charged particle polarizabilities can be 
done directly on the lattice by measuring a Compton scattering coefficient [To| ]. However, the direct evaluation is 
actually quite involved because of the many disconnected diagrams involved. There is a way around this difficulty if 
one is willing to work in the exact chiral limit by use of the Das, Mathur Okubo (DMO) sum rule Although this 
results in a simpler lattice calculation, a well-known problem in this approach is that the spectral integral involves a 
rather precise cancellation of large numbers and so may be difficult to evaluate numerically. This paper attempts to 
understand to what extent this difficulty holds in lattice evaluations and to begin to explore the internal systematics 
of such calculations. 

An explanation of the methods used to estimate the DMO spectral integral will be given in the next section. 
Successful use of the quark continuum model will be made; this sheds additional light on it's use and range of validity. 
The consistency of the results will be probed using an alternate purely numerical approach based on Simpson's rule. 
We will also check the consistency of our data with low energy phenomenology by extracting the pion and rho meson 
decay constants, / w and / p , from the axial and vector correlation functions. Various other sources of systematic error 
will be estimated and suggestions for further improvements in future lattice evaluations of charged pion polarizability 
will be made. 

II. Background 
A. Review 



The DMO sum rule for the charged pion polarizability, derived from current algebra in the chiral limit is: 

= a^2L - -2L_ j°° ^W 2 ) - PA{s 2 )). (1) 



3m w 2m T /2 



(m-Tr+e) 2 



The pion pole is absent from the axial vector spectral integral; this is symbolized by the lower limit, (m w + e) 2 , e being 
a small positive quantity. A conventional lattice normalization for the spectral functions pv,A is being used here; see 
Eq.(3) below. 

There are several advantages to using Eq. (|l]) instead of attempting direct evaluations of pion polarizability. Eq. (Q) 
employs only mesonic two-point functions, which are easily calculable on the lattice. In addition, the necessary 
propagators can be precisely measured using a combination of numerical techniques which will be explained in Section 
IIIA. 

There is no unique way of evaluating Eq.(|l|) on the lattice. The strong interaction parameters, pion mass (m^), 
pion decay constant (f n ) and the charged pion squared charge radius (< r\ >) also enter this expression. All these 
quantities can and should have separate lattice evaluations. In the limited point of view adopted here, evaluation of 
the spectral integral in Eq.(|l|), the so-called intrinsic part, 

alt ^ r ^-( Pv ( S i) - p A {s*)l (2) 



will be concentrated on. We will also compare the implied values for the pion decay constant f v and vector meson 
decay constant f p against experimental results. Normally since polarizability scales like a 3 , one expects it to be very 
sensitive to the lattice scale. However, because experimental values of and /„• will be used in the DMO expression 
for a™t there will actually be no scale uncertainty in the answer obtained here. The only quantity for which we will 
need the lattice scale will be in the later calculation of /„■. 
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B. Derivation 

The Minkowski space DMO sum rule spectral densities pv(s 2 ) and pa{s 2 ) arise from jl4| 

T[v 3 o (x)v 3 o H0)~al(x)a^(0)}\0) (3) 
1 f ds 2 2 \ / 2\\ ( kfj,k v \ f 2 k^k v 
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(p y (s ) - p A (s )) -g^ + 



(m„+e)= fc 2 + S 2 V s / fc 2 +/i 



where fco = \/s 2 + fc 2 or fco = y /i 2 + fc 2 (/i= pion mass) in the pion term. Also in Eq.(3) the metric is = 
(+, — , — , — ), v^(x) = -0(a;) ; ^-7 ^ip(x) , and a°(x) = ^{x) 1 ^-^^^^). This same amplitude can also be expressed as 

i J d^xe-^Ol T K 3 (x) Uo 3t (0) - al(x)a^ (0)]|0) (4) 

=iy ^(oiTK d (^< d t (0) _ a ^ (a;)a ^t (0)] | 0) . 

in the exact 5C/ (2)-fl.avor limit because of either canceling or vanishing self-contractions of the currents. The currents 
being used in Eq.(4) are v^ d {x,t) = d(x,t)jou(x,t) and a,Q d (x, t) = d(x,t)"f5"fou(x,t). It is noticed that 

( pv ( s *)-p A ( s *))=i / dKA&(* 2 ,t)-A&(k 2 ,t))| Sa=0 , (5) 

where 

A£,(fc 2 ,t) = |d 3 x e -^(0|TK d (x,tK dt (0)]|0), (6) 

A^fc 2 ,*) = J d 3 xe-^(0\T[ar(x,t)a u o d H0)}\0). (7) 
Last of all, a switch is made to imaginary time, 

dtA 0Q (k 2 ,t) = dt E A A4 (k 2 ,t E ), (8) 



where we have the standard ud-flavor propagators (real and positive for t ^ 0) 

A^(k 2 ,t E ) = J d 3 xe~^ s (0\T{vf(x,~it E )vf\0)}\0), (9) 
Af 4 (fc 2 ,^) = [ d 3 xe- lts (0\T[af (x,-it E )af\0)}\0). (10) 



(One may take 70 = 74 here.) Thus, putting Eqs.(|2j)-([10|) together, we have 

<± = -^Tf2 3^ r dt E (AUk 2 ,t E )-At 4 (k 2 ,t E ))\^ 2=0 . (11) 

Zm TTj7T dk Z J -OO 



On the lattice, the right-hand side of Eq.(|ll|) will be formed by taking a numerical momentum derivative of the 
lattice propagator data. Although the derivative in (|^) removes the pion contribution in Eq.(3), one has only finite 
momentuma on the lattice and the pion contribution must be explicitly subtracted from the axial vector propagators 
in the lattice version of Eq. dll|) . 

C. Continuum Model 



The time integral in Eq. (11) will be performed in both a discrete and continuous sense using the lattice data; the 
difference will be taken as a reasonable estimate of the systematic error of the integral. In the discrete case, one 
can simply apply Simpson's integration rule to the Euclidean time propagator data after subtracting out the pion 
contributions. This will be described in detail later. The continuous case demands some way of interpolating between 
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the propagator time values. For this purpose, let us consider the standard lattice Euclidean ud-flavor point-to-point 
charge density correlator (<f = ka): 



J2 e~^ S < 0\T(vf(x)vf\0))\0 > . 

X 

One may show that this reduces to (replacing the sum ^Tg by the integral J d 3 x) 

d 3 x e- l **Tr[S(x, 0) l4l5 S\x, 0)7475], 



where S(x,y) is the quark propagator and the trace is over color and Dirac spaces. Using the free quark propagator 
(coordinate gauge; diagonal to this order in color space), 

S( a;> 0) = -Ll^ + 7 -^^ + ..., (12) 

the following definition (q = \q\) is now made: 

G 44 (t,q) = J d 3 xe-^ s £- r t * r (13) 

Doing the angular integrals yields 

12 f°° 2t 2 1 

G 44 (t,q) = — drr$m{qr){ - }. (14) 

Actually, what one wants in this case is a derivative of the above with respect to squared spatial momentum (see 
Eq. ( pi] ) above) evaluated at zero momentum. Since only finite momentuma are available in the lattice simulation, this 
continuum procedure can not be reproduced on the lattice. Thus in order to compare with the lattice data consider 
instead, 

A ^ G 44 (t,q)-G 44 (t,0) ^ 

Q 



where q represents the lowest lattice momentum value. Fitting the lattice data with Eq.([l^) has the advantage of 
using the same type of "derivative" as in the lattice data, but has the disadvantage of including a small momentum 
dependence in the phenomenological model. To this order the same functional form holds for the axial vector propa- 
gator as well. The next nonvanishing term in the lattice propagator from Eq.jl^) is proportional to the square of the 
quark mass, but this term gives negligible contribution to the fits and is not considered further. 

A number of modifications are necessary to this function before one compares to the lattice data. First, it is clear 
that the above expressions have an ultraviolet infinity associated with the r — lower limit for t = 0. This infinity 
can be controlled, as the lattice itself controls it, by putting in a short distance cutoff. So, replace the lower limit 
above by ro > 0, which becomes a parameter in the fits. Call this modified function AG 44 (t, q, Tq). Second, put in a 
continuum threshold, So, to control the onset of excited states in the spectral density. One can show that the resulting 
function is given by 

/•oo roc 1 

Ag 44 (t,q,r ,s )= ds —e ms AG 44 (\t\+iu,q,r ). (16) 

J s J -oo 27T 

A third parameter, £, will be introduced as a multiplicative factor normalizing this function. Originally, this parameter 
was introduced to account for lattice anisotrophy at small lattice time separations (see Ref. pH). However, for an 
interesting alternate interpretation of this parameter, see Ref. [0. One can analytically perform the two integrals 



in Eq.(16); the remaining oscillating radial integral is done numerically. This numerical evaluation makes the fitting 
of parameters for this model rather slow and renders a third-order jackknife analysis of the data impractical. For a 
more explicit representation of the function defined in Eq.(|l^), see the Appendix. 

In attempting to fit the (momentum differentiated) propagator time data, pole mass terms of the form 

(X p fe~ E ^ + {X v ) 2 e- Evt , (17) 
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will be added in the vector case and the three terms 



(XI fe- E - 1 - (A°) 2 e- m "' + {\ q A fe- EA \ (18) 

will be added in the axial case. The E p term (in (|l7j)) and the E„, terms (in (Q)) represent the lowest pole 
contributions to the spectral integral in the vector and axial vector channels. (There is no term proportional to e _, ™ p * 
in the vector case because of charge conservation.) The additional pole terms in these expressions, the Ey term in 
( |l7j ) and the Ea term in (18), were found to be necessary in achieving a reasonable fit to the lattice data; see the 



comments in Section IIIB. The mass values and m p will be fixed from independent lattice measurements, and the 
continuum relationship 



E^a = ^q 2 + (m,a) 2 . (19) 

will be assumed. Phcnomenologically, the Ey pole in the vector case has quantum numbers of radially excited states 
the rho meson and the Ea pole has the quantum numbers of the a\ meson. All told, there are 6 parameters in the 
vector fits (£y, soy, roy,X q p , X v , Ey) and 7 parameters (£4, sq.a, tq^a, A°, A|, X q A , Ea) in the axial fits. (However, 
see later comments about the continuum relation between A° and X% in Section IIIC.) Of course, once the axial data 
is fit, one must eliminate the pole terms involving the pion before doing the DMO integral. 

In the context of the continuum plus pole model, the Euclidean lattice data is now fit to the forms, 

A^ 4 (q 2 , t) - A&(0, t) = q 2 £ v Ag 44 (t, q, r , s ) + [Xpfe^ + [\y) 2 e- Ev \ (20) 



t) - A&(0,t) = q 2 UAg^(t,q,r ,s ) + (X^e 



\2„-m 7r t 



-E A t 



(21) 



A 44 (q 2 , t) and A 44 (q 2 ,t) represent the actual correlation functions measured on the lattice. It is assumed that these 
are related to the continuum functions needed in Eq.([ll|) by scale, tadpole and renormalization factors as follows: 



a 3 AX 4 (k 2 ,t) 



jv7 



(1 -0.82a y ) 



A&(g 2 ,t), 



(22) 



a 3 At 4 (k 2 ,t) 



~(l-Q.31a v y 



(23) 



where ay is the strong interaction coupling constant defined in Ref . |Tqj , N s is the number of spatial sites smeared 
over in the source interpolation fields (16 2 in this case), and 



Nt = 16(1 — 



(24) 



We will use cty(— ) = 0.1557 (j3 = 6.0) for the local currents in this study. In the following tables and figures, the 
results for the intrinsic part of the charged pion polarizability will presented in natural dimensionless form, Ty t A, 
where 



1v,A z - 



Q 



m^f 2 



(J A — T\, 



ds 2 



Pv,a{s 2 ). 



(25) 
(26) 



III. Results 
A. Simulation Parameters 

The simulation was done on 32 quenched configurations with Wilson fermions on a 16 3 x 24 lattice at (3 = 6.0. 
The lattices were constructed with the algorithm of Ref. thermalized by 11000 sweeps and separated by 1000 
sweeps. Four values of the Wilson hopping parameter were considered, n = 0.154,0.152,0.150 and 0.148. We used 
the "Volume method" pl[ to calculate the propagators directly from the nongaugefixed configurations, smearing over 
a 16 x 16 spatial plane at time step 8 of the lattice. (The first time step of the lattice will be defined to be t = 1.) 
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This sacrifices the Fourier transforms in two spatial directions but reinforces the momentum projection in the third. 
This is the main idea of the "Fourier reinforcement" method [E^] . 

The signals obtained for the vector and axial vector charge density operators are excellent. Figs. 1 and 2 show the 
local lattice energy, 

for the vector and axial vector cases, respectively. The p and n masses were calculated separately with extremely long 
time baselines on 20 of these configurations, fixed to the lattice Coulomb gauge, with quark propagators starting at 
t = 1 and single exponential fits to time steps 16 to 19. The results are given as the first line in Tables II and III. 



Three of these results (k = 0.154, 0.152, 0.148) are taken from Table I of Ref. |23|; the result at k = 0.150 is new. The 
value of K cr = 0.1564 is also taken from Ref. (^3|. The dimensionless quark mass in this reference as well as here is 
taken to be 

ma = Z n (^_ 3 ). (28) 

The correlated chi-squared per degree of freedom, x% on all of the mass fits (no SVD decomposition; see Section IIIB), 
were less than one. The horizontal lines in Figs. 1 and 2 show the predicted energies using continuum dispersion; 
agreement with the lattice data on time steps 15 to 20 inclusive is evident. These masses will be used as fixed input 
rather than parameters in the present calculation, which significantly improves the error bars on the remaining fit 
parameters. The systematic effects of varying the input p and 7r masses in this calculation will be reported on in 
Section IIIC. 



B. Spectral Integral Evaluation 

The data in the vector and axial vector sectors was fit to Eqs.(^0|) and (pl|); examples of these fits are shown in Figs. 
3 and 4, which show the case n = 0.154. It is the area under these curve which is of interest here. The time integrals 
are very sharply peaked and the fits themselves extend from the source timesite at t = 8 to t = 20. The numerical 
quality of the vector time data is seen to be better than the axial vector case, but both are quite acceptable. The 
parameters of these fits are given in Tables II and III. 

A second-order single elimination jackknife was used for error analysis at each k value; the first order defines error 
bars on the time correlation functions and the second defines errors on the fit parameters of these functions. The 
fits reported in Tables II and III are characterized by their correlated chi-squared per degree of freedom, x% These 
were arrived at by the singular value decomposition (SVD) algorithm suggested in Ref. p4[ . It was found there that 
the correct correlated x\ was obtained on small data samples when the number of exact eigenvalues retained, E, was 
chosen to be sw y/N, where N is the number of configurations. We will use E = 6 (N = 32). Defining D to be the 
number of fit timesites (D = 13), Ref. finds the increase in x\ from the "true"result is given by 1 + (D + 1)/N, 
which in the present case is approximately 1.44. (Comparing SVD and non-SVD fits, this ratio was actually found to 
be 1.57 in the vector case and 1.45 in the axial case, averaged over k.) The x\ values in the vector case (~ .3 — .6) are 
quite good; the x? values in the axial case (~ 1.0— 1.5) were higher, but are still acceptable. The pole terms involving 
Ay, Xa in Eqs . (pp|) and ( p0| ) were crucial to obtaining acceptable fits in both the vector and axial vector sectors. The 
fitting of 13 timesites with acceptable x\ values including the time origin is an extremely nontrivial matter and shows 
the usefulness of the (cutoff) continuum quark model in fitting lattice propagator data. 

As is suggested by the numerical results of Ref. pij , the best values of the fit parameters (and their error bars 
from the jackknife) are actually determined by doing uncorrelated fits; the correlated x\ was use d only as a selection 
criterion of fit time intervals. This is the same procedure as used in Ref. p^| . The error bars in the parameters of 
the time fits were determined by the jackknife, while the errors in the fits across k values were determined by the 
Levenberg-Marquardt method |2C| using the CURFIT routine of Ref. (27]]. 

Tables IV and V give the relative contribution of the various continuum model sectors to the final result. One sees 
that the (subtracted) axial sector is almost saturated by the assumed pole, but that only about 30% of the vector 
result is given by the p-meson pole, the majority coming from the continuum. This is very different from the chiral 
model expectations and has important consequences for the final answer. 

The fit parameter values in Tables II and III are fairly reasonable. The continuum threshold values So.Vj $o,A are 
of order unity, with so,y tending to cluster just above the p-meson mass; there is no particular trend in the axial case. 
(Large fluctuations in so,y occur at k = 0.154 in the vector case although the jackknife error in the final integral 
remains quite small; see Table VII.) The Ey, Ea pole energies, mimicking the contribution of higher bound states, are 
seen to decrease with increasing k. The values of A p , A^ and A* are also reasonable and will be examined extensively 
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in Section IIIC. However, there are also some questionable aspects to the parameter values in Tables II and III. In the 
interpretation of Ref. fig] , the £,v,A values should be approximately constant across k. In fact, there is an increase in 
these values as k increases in these tables, although the axial case is dominated by errors. (Note that the £,v,A values 
are significantly decreased by the inclusion of the Xy,A pole terms.) In addition, the short distance cutoff values r^y 
and tq^a are not particularly constant in k as one might expect. The axial case Tq^a values are dominated by errors 
and no real comparison between the two sets of values can be made. 

As a completely independent means of approaching these integrals, the Simpson integration formula for discrete 
data [p6| was utilized. In order to do this, it is necessary to explicitly remove the pion poles in the zero momentum and 
nonzero momentum axial vector propagators. These were fit with single exponentials (using the SVD decomposition) 
to the 18 to 20 time sites, with acceptable correlated \\ values in both cases. Table VI reports the results of these 
evaluations. The A°, \% values reported in this Table are then used to remove the pion tails from the axial vector 
propagator time integrals. There is no unique way of numerically integrating the subtracted data. Since the lattice 
data is so strongly peaked in time, different integration rules can give significantly different results. The Simpson 
rule was chosen because it was simple and well-known, but other rules could have served as well. The choice of an 
odd number of time sites to fit (13) means that the discontinuities in the fit polynomials will occur at odd time sites, 
including the origin. This allows for the strong peaking near the time origin seen in Figs. 3 and 4 and gives much 
better agreement with the continuum model integral results than an integral rule which requires continuity at the 
origin. 

Fig. 5 shows the final chiral extrapolations in k, with numerical results reported in Table VII. This Table shows 
that the Simpson and continuum results for Ty are rather close, whereas 1a is the major source of the systematic 
error in the time integral. The statistical errors in I a are also significantly larger than those in Ty. The continuum 
model values are shown as squares, the Simpson evaluations as circles, with the filled-in symbols representing the 
extrapolated values. The Simpson values are larger but extrapolate to a smaller result because of the positive slope. 
On the other hand the continuum model results have a more sedate, negative slope. Combining these results, the 
dimensionless integral X is now given as 

J= 36.3(3.9)(3.5) x 1CT 3 , (29) 
implying (using = 139.6 MeV, f„ = 92.4 MeV) 

a rnt = „ 17 . 1 ( L8 )( 1 . 6 ) x 1(T 4 /™ 3 , (30) 

where the first number in parentheses is the statistical error, the second is the systematic error, taken to be half of 
the difference of the continuum and Simpson model central values. The central value and the statistical error are the 
average of the continuum and Simpson results from Table VII. The central value in Eq.(|29|) is significantly larger than 
current experimental or chiral results would imply. For example, Ref. [[| quotes 27.0(.5) x 10~ 3 and Ref. ||] gives 
21.0(.5) x 10~ 3 for this same quantity. When combined with the experimental result < r\ >= 0.439(.008)/m 2 from Ref. 
||§, this implies a negative pion polarizability, a I± = -2.0(1. 8)(1.6) x 10~ 4 /m 3 . (The result < r 2 . >= 0.463(.006)/m 2 
from Ref. |2!| would imply a,± = — 1.2(1.8)(1.6) x 10 _4 /m 3 .) We will explore the systematics responsible for this 
outcome in the next subsection. 

C. Systematics 

The pion and rho meson masses listed as the top lines in Tables II and III are themselves measured from the lattice 
and have their own Monte Carlo statistical errors. Since these are treated as input rather than fit parameters, one 
should investigate the systematic errors associated with varying these inputs. This is done by studying the change in 
the central value of the chiral-extrapolated result, Eq. fl29|) , when the rho and pion masses are put at the upper and 
lower limits in Table III. When this is done with the rho meson, the central value changed by approximately ±1%. As 
one might expect, the Eq. (p9|) result is more sensitive to the input pion mass because this is used to fit and remove 
the pion spectral contribution. It is found that the results from the continuum quark model are considerably more 
sensitive to the input pion mass than the Simpson results. By varying the pion mass within the limits in Table III 
one finds that the central value of the result in Eq.(|2^) changes by (symmetrizing the upper and lower changes) by 
about ±2.3%. This is small compared to the estimated systematic uncertainty due to the integral model (continuum 
or Simpson) dependence. 

In order to estimate the size of the finite lattice spacing errors in this simulation, we replaced the continuum relation 
Eq.(|l9|) with the lattice spin dispersion relation |Q) 

sinh\^) = sinh^) + £ sm 2 (^), (31) 
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and made the substitution |3lJ], 

g — 2«in(|), (32) 

everywhere for the momentum factor. When these changes were made, the central value in Eq.(|2^) was increased by 
approximately 2.5%. 

As a consistency check of the correlation functions used in this study with low energy phenomenology, the values 
of the vector meson and pion decay constants have be measured and compared to experiment. The vector meson case 
will be especially revealing since this is a dimensionless quantity independent of the lattice scale. 

The continuum matrix element for rho meson decay is given by 

(0|<(0)|p A (p)) = ^^= £ >,A), (33) 



where f p is the decay constant, the polarization state is specified by A and /i = for the charge density operator. One 
can show that in terms of the A p parameter in the continuum fits, the implied vector meson decay constant is given 
by 

U = q{mVa) ■ (34) 

where q = ~. The results for the vector meson decay constant, f p , are given in Table VIII and are shown in Fig. 6. 
When all 4 k values are extrapolated to K cr , we obtain /„ = 3.11(.39) (Xd = -59), lower than the experimental value 
of 3.56(.14). However, it is well known in lattice studies on similar sized lattices that the p to nucleon mass ratio, 
mp/mpj, is underestimated when extrapolated to the chiral limit. Using the nucleon mass to set the lattice scale, 
if one instead extrapolates f p to the p-meson physical mass, about 770 MeV (occurring at about k = 0.1545), one 
obtains a much better result, f p = 3.48(.28). We will check that the f n values from the axial vector propagator given 
by this choice of scale is consistent with experimental results. 

The large time limit of the axial vector propagators imply values of the pion decay constant, /„•. One has that 

(0|<(0)Hp))=*^, (35) 



where again p, = 0. One can relate the parameter A° to the pion decay constant as follows: 

U = — (36) 

The results of this calculation are shown in Fig. 7 {or 1 = 1.74 GeV). The extrapolation to the chiral limit is extremely 
straight and one obtains 87.4(9.0) MeV, consistent with the experimental result of 92.4 MeV. 

As a check on the /„■ calculation, we consider the ratio of the axial pole parameters for the pion, which by the 
continuum relation Eq.(pq) is given by 



A«/A° = VE~/nw. (37) 

This comparison is carried out in Table IX for the continuum and Simpson fits. There appears to be a small violation 
of this continuum relation at perhaps the 5 — 10% level in the lattice data. However, note that a systematic error in 
this quantity does not necessarily affect the calculation since the pion contributions are excluded from the DMO axial 
integral. Also note that the two data treatments (continuum and Simpson) are quite consistent with one another for 
this ratio. 

We have seen above that f p , extrapolated to the chiral limit, gives a result which is small compared to experiment but 
that the / w value extracted using the axial propagator and the nucleon mass scale is consistent with phenomenology. 
This suggests that it may be the behavior of the vector propagator in the chiral limit which is responsible for the 
large value of X in this calculation. To test this, one may instead extract the value of the vector contribution Xy at 
the physical p-meson mass, similar to what was done above for f p ; we now obtain X = 33.2(4.0)(2.9), a better result 
but one which is still too large. The remaining difference is clearly due to the large "continuum" contribution to the 
vector correlator remarked on earlier; see Section IIIB and Table IV. 

Finally, note that the intrinsic polarizability, through the renormalization factors in Eqs.(p^) and (|2^), is fairly 
insensitive to the value of ay used since the vector and axial vector integrals contribute with opposite signs. Roughly 
speaking, it is found that a x% percentage change in the value of ay induces a change in Eq.dga) of about —x/2%. 
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IV. Conclusions and Remarks 



The sum rule method of extracting charged pion polarizability from lattice data has been examined. In the limited 
point of view adopted here, the spectral integral in the DMO sum rule has been considered separately. It has been 
evaluated with lattice data and the result, large compared to phenomenology and chiral models, is given by Eq. (^9|) . 
Smaller systematic effects from altered input mass values, finite lattice spacing and renormalization constants were 
also considered. 

Excellent fits were obtained to the (momentum-differentiated) lattice data across 13 time slices, including the 
propagator origins, using the continuum quark model. The time fitting of these quantities with reasonable \\ values 
would not have been possible without: 1) introduction of the lattice cutoff, fp; 2) addition of additional pole terms in 
both the vector and axial vector cases; 3) the SVD modification, following Ref. |2q] , of the time propagator eigenvalues. 
We have also seen that the axial vector propagator is largely responsible for both the statistical and time-integral 
systematic errors. 

The lattice systematics have been examined extensively and it has been argued above that 1) systematics associated 
with the incorrect lattice ratio m p /mN and 2) the fact that the vector propagator is far from being dominated by 
the p- meson are responsible for the large central value of our final result in Eq. (^6|) . It is precisely because pion 
polarizability is so sensitive to and revealing of lattice systematics that it represents a significant test of the ability 
of the lattice to produce phenomenologically interesting predictions. Further studies with larger lattices and better 
actions should be even more revealing of these systematics. 

It is clear that in order to obtain phenomenologically interesting values of a 1 ™* from the lattice, both the statistical 
and systematic errors here will have to be reduced. The systematic uncertainty in the time integrals can be reduced 
by using a time-asymmetric lattice with a fine mesh of lattice points in the time direction. This will allow sampling 
at smaller time intervals (but not too small to get into the asymptotic time regime) in evaluating the strongly peaked 
integrals. However, in order to understand the dynamics leading to pion polarizability it will be necessary to go beyond 
the DMO sum rule to direct measurements. These additional considerations will be taken up in future publications. 
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Appendix 

In this brief Appendix, a more explicit form of the continuum model expression for AQ^t, q, ro, so) will be given. 

In introducing the ultraviolet cutoff, one can make use of the procedure in Ref. Jl6| , which gives an upper limit, A, 
in energy. The resulting expressions are rather complicated and are only defined in the limit t — > 0. Alternatively, 
one can define a short distance cutoff simply by putting a lower limit on the r-space integrals. This is the procedure 
followed here. Then combining Eqs.(^) and ( |l5| ) of the text, we have 

A _ , , 12 f°=, ,sin(qr) 2t 2 1 

AG44M)ro) = __^ drr( ^__ r){ _^ ____}. (38 ) 

The continuum threshold that is introduced in Ref. []l6f is equivalent to doing an an incomplete Laplace transform of 
the spectral density, which is itself obtained with an inverse Laplace transform of the propagator. A similar procedure 
is followed here with the vector, axial vector continuum spectral densities. The basic assumption of the continuum 
model is that we are at low enough lattice momentum that 

2 n r°° rlv 2 

AG 44 (t,g,ro)^° J gfe-'M* 2 ) (39) 
is a reasonable identification. In this limit this means that AG44(£, ro) has the time integral 

/OO 2 /"OO j 2 

dtAG M (t 7 q,r ) q ^° ^-P(s 2 ), (40) 
-00 Jo s 

consistent with Eq.(3). Now putting Eq.([38|) into Eq.(^6[) results in the explicit expression (for t > 0), 
/ n 12 r ^ /-°° , ,sin(qr) . [°° , f°° , e ms (t + iu) 2 



dr r( S ^ n gr — r) / c?s / du— . , 2 x 2 }- 



9 •/-«, (r 2 + (t + m) 2 ) 5 

The poles in the M-integral are identified and the integral done. The remaining s-integral is then done explicitly. The 
final result can be presented as follows: 

AG^q,r , So) = ' Re{ f -(^^ - r)e^ ){ « + (42) 
in q J ro r q [t — iry (t — ir)° 

3sq si 3.1 so , 3i 
H 1 1 I 1 1 H 11 

(t - ir) 2 (t- ir) r 2l (t- ir) 2 (t - it) 1 r 3 (t- ir) U 

3 r r t dr f sin(qr) _\ o a (— *+*»■) r 

2 2s 



Im{ / y"v»-' _ r )e s ^- t+4r, { 



2ir 3 q 2 J ro r 2 q (t — ir) 3 (t — ir) 2 

^--{7^ + 77^} + !^}}- 



(t — ir) r (t — ir) 2 (t — ir) r 2 (t — ir) 

The explicit real and imaginary parts are then separated out from this expression and the remaining r-integral is done 
numerically in the continuum fits of the vector and axial vector data. 
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TABLE I. Previous experimental and theoretical results on charged pion polarizability. 



a n ± (units: 10 4 fm 3 ) Type Reference 



6.8 ± 1.4 ± 1.2 Expt. Ref. [1 

20 ±12 ±1.2 Expt. Ref. J \ 

2.2 ±1.6 ±1.2 Expt. Ref. |j 

2. 64 ±.36 Theor. /Expt. Ref. [5 

~ 3.6 Theor. Ref. \\ 

2.4 ±0.5 Theor. Ref. [7 

5.6 ± 0.5 Theor. Ref. | 



TABLE II. Continuum model vector fit parameters. \d gives the chi-squared per degree of freedom for the fit. 



Quantity 


k = 0.154 


0.152 


0.150 


0.148 


m l p nput a 


0.463(.020) 


0.550(.013) 


0.635(.010) 


0.718(.008) 




4.44(.54) 


4.16(.78) 


3.71(.52) 


3.57(.50) 


so,v 


0.90(.55) 


0.568(.050) 


0.644(.050) 


0.724(.049) 


ro,v 


0.442(.028) 


0.417(.040) 


0.383(.029) 


0.368(.029) 




4.86(.43) x 10" 2 


4.36(.41) x 10" 2 


4.65(.28) x 10" 2 


4.79(.18) x 10 


\v 


4.6(1.7) x 10" 2 


4.6(1.1) x 10" 2 


5.03(.62) x 10" 2 


5.12(.60) x 10 


Eva 


1.17(.14) 


1.33(.18) 


1.46(.10) 


1.55(.07) 


x\ 


0.34 


0.30 


0.44 


0.60 
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TABLE III. 


Continuum model axial fit parameters. Xd gives the chi-squared per degree of freedom for the fit. 


Quantity 


k = 0.154 


0.152 


0.150 


0.148 


m*" put a 


0.366(.010) 


0.479(.008) 


0.581(.007) 


0.676(.006) 


U 


1.5(1.7) 


1.1(1.2) 


0.96(.45) 


0.61(1.0) 


SO, A 


1.19(.24) 


1.31(.18) 


1.34(.17) 


0.69(.23) 


ro,A 


0.50(.30) 


0.49(.24) 


0.546(.082) 


0.60(.54) 


\0 

7T 


3.73(.26) x 10~ 2 


4.96(.25) x 10~ 2 


6.18(.26) x 10~ 2 


7.37(.25) x 10~ 2 


\% 


5.05(.37) x 10~ 2 


6.10(.34) x 10~ 2 


7.16(.34) x 10~ 2 


8.15(.41) x 10~ 2 


A A 


7.6(1.3) x 10~ 2 


7.95(.84) x 10~ 2 


8.04(.38) x 10~ 2 


8.22(.85) x 10~ 2 


Eaci 


1.38(.15) 


1.49(.ll) 


1.58(.06) 


1.70(.07) 


Xd 


l.UU 


l.Oo 


1.40 


1 1 Q 




TABLE IV. 


Continuum model vector 


integral strengths. 




Contribution 


k = 0.154 


0.152 


0.150 


0.148 


Continuum 


0.51(.18) 


0.61(.09) 


0.56(.06) 


0.55(.05) 


Lowest pole 


0.34(.06) 


0.25(.05) 


0.27(.03) 


0.28(.02) 


Excited pole 


n i c\( i 
U.lO^.lo ) 


0.14(.05) 


0.16(.03) 






TABLE V. 


Continuum model axial integral strengths. 




Contribution 


k = 0.154 


0.152 


0.150 


0.148 


Continuum 


0.28(.25) 


0.20(.17) 


0.18(.08) 


0.16(.25) 


Pole 


0.72(.25) 


0.80(.17) 


0.82(.08) 


0.84(.25) 
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TABLE VI. Simpson model parameters. \\ gives the chi-squared per degree of freedom for each fit. 



Quantity 


k = 0.154 


0.152 


0.150 


0.148 


A" 


3.70(.24) x 1Q-' A 


4.94(.22) x 10~^ 


6.23(.23) x W~~ z 


7.53(.24) x 1Q-~ Z 


xl 


0.57 


0.83 


0.82 


0.80 




4.73(.31) x 10~' z 


5.95(.27) x 10~^ 


7.27(.27) x 1Q~ Z 


8.63(.28) x 1Q-~ Z 


Xd 


0.19 


0.15 


0.13 


0.09 



TABLE VII. Results for the dimensionless integral X; a factor of 10 3 multiplies all the entries. The "C" superscript indicates 
the continuum model values and the superscript "Syndicates values from Simpson fits. 



Quantity 


K cr = 0.1564 


k = 0.154 


k = 0.152 


k = 0.150 


k = 0.148 


7 C 

L v 


79.1(1.9) 


74.6(1.5) 


72.1(1.3) 


68.5(1.3) 


65.1(1.2) 


7 s 
L v 


79.8(1.2) 


76.4(.9) 


74.8(.9) 


72.3(.8) 


69.4(.7) 




39.6(4.8) 


37.1(4.0) 


34.6(2.6) 


32.4(1.5) 


30.7(3.6) 


-La 


47.0(2.9) 


40.7(2.5) 


36.8(1.9) 


32.1(1.8) 


27.1(1.7) 


l u 


39.8(4.6) 


37.5(3.9) 


37.6(2.5) 


36.1(1.6) 


34.3(3.5) 


X s 


32.8(3.1) 


35.6(2.7) 


37.9(2.0) 


40.1(1.9) 


42.3(1.8) 



TABLE VIII. Vector and pion decay constants. f£ indicates values inferred from continuum model fits and f% values from 
the Simpson fits. 



Quantity 


k = 0.154 


0.152 


0.150 


0.148 


f P 


3.40(.30) 


4.26(.40) 


4.39(.26) 


4.60(.17) 


f° (MeV) 


107(7) 


125(6) 


141(6) 


156(5) 


(MeV) 


106(7) 


124(6) 


142(5) 


159(5) 
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TABLE IX. Comparison of the continuum ratio AJ/A° for continuum and Simpson model axial fits. 



Quantity 


k = 0.154 


0.152 


0.150 


0.148 


Continuum theory 


1.21 


1.14 


1.10 


1.08 


Continuum model 


1.35(.10) 


1.23(.06) 


1.16(.05) 


1.11(.04) 


Simpson model 


1.28(.09) 


1.20(.06) 


1.17(.04) 


1.15(.04) 
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Figure Captions 



1. Local E l p att (t + i) measurements for the lattice vector charge density, A v (q 2 ,t), versus lattice time location, 
t, compared with continuum dispersion (horizontal lines) for |<f| = |. Squares are for n = 0.154, triangles are 
k = 0.150, stars are for k = 0.148, and circles are k = 0.148 results. 

2. Local E l ® u (t+7j) measurements for the lattice axial vector charge density, A A (q 2 1 t), versus lattice time location, 
t, compared with continuum dispersion (horizontal lines) when \q\ = |. Symbols are the same as in Fig. 1. 

3. The base 10 logarithm of the difference in nonzero and zero momentum lattice vector charge densities, 
logio(A v (q 2 , t) — A y (0, t)), at k = 0.154 as a function of lattice time location. The continuum model fit 
is shown. 

4. The base 10 logarithm of the absolute value of the difference between nonzero and zero momentum lattice axial 
vector charge densities, logio\A A (q 2 , t) — A A (0,t)\ at k = 0.154 as a function of lattice time location. The 
continuum model fit is shown. Notice the change in sign of the correlation function slightly after time location 
12. 

5. Chiral extrapolation of the final results for the dimcnsionlcss spectral integral, X, as a function of dimcnsionlcss 
quark mass, ma. Squares are for the continuum model results and circles are for Simpson model results from 
Table VII. The filled- in symbols represent the chiral-extrapolated results at ma = 0. Note that some of the 
circles are offset in ma for clarity of presentation. 

6. The vector meson decay constant, f p , as a function of dimcnsionlcss quark mass, ma. The boxes indicate 
measured values at k = 0.154,0.152,0.150 and 0.148. The filled circles indicate the extrapolated result at 
ma — as well as the extrapolated result at the physical p meson mass. The dotted line shows the experimental 
upper and lower limits on f p . 

7. The pion decay constant, /„., as a function of dimcnsionlcss quark mass, ma, in MeV. The boxes indicate 
measured values at k = 0.154,0.152,0.150 and 0.148. The filled circle indicates the extrapolated result at 
ma = and the dotted line is the experimental result, f v = 92.4 MeV. 
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